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Abstract 

Using diagrammatic methods, we show how the Ward identity can be used to constrain the ladder 
kernel in transport coefficient calculations. More specifically, we use the Ward identity to determine 
the necessary diagrams that must be resummed using an integral equation. One of our main results 



. is an equation relating the kernel of the integral equation with functional derivatives of the full 

> . 

I/-) , self-energy; it is similar to what is obtained with 2PI effective action methods. However, since we 

CN ' use the Ward identity as our starting point, gauge invariance is preserved. Using power counting 



arguments, we also show which self-energies must be included in the resummation at leading order, 
including 2 to 2 scatterings and 1 to 2 collinear scatterings with the LPM effect. We show that 
Q_i| our quantum field theory result is equivalent to the one of Arnold, Moore and Yaffe obtained using 



effective kinetic theory. In this paper we restrict our discussion to electrical conductivity in hot 
QED, but our method can in principle be generalized to other transport coefficients and other 
theories. 
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I. INTRODUCTION 



Transport coefficients are measures of the efficiency at which a conserved quantity is 
transported on "long" distances (compared to microscopic relaxation scales) in a medium. 
For example, electrical conductivity characterizes the diffusion of charge due to an external 
electric field and shear viscosity characterizes the diffusion of momentum transverse to the 
direction of propagation. In real non-relativistic systems (e.g. interacting electron gas in 
a lattice), transport coefficients are almost impossible to calculate, because systems are 
strongly interacting and no simple closure exists; but in hot, weakly interacting theories, 
they can in principle be evaluated. The computation of these quantities is important from 
a theoretical point of view: since they characterize linear deviations away from equilibrium 
but are computed from well-known equilibrium field theory tools, they could be used as 
a benchmark for testing non-equilibrium field theories, which are less well developed than 
equilibrium ones (but see the recent developments in Refs. They could also 

have an influence^pn^ the physics of the early universe, such as electroweak baryogenesis 
(see for examp 



example 



6 



7 



e {4, 5|) and the formation and decay of primordial magnetic fields (see for 
). Shear viscosity also attracted a lot of attention lately in the heavy ion 



community, partly due to the results on elliptic flow (seemingly implying a low viscosity, 
see for example p, []l3] and the references therein) and its exact computation in a strongly 
coupled Super Yang-Mills theory 11] . The above examples show the importance of having 



a good theoretical handle on transport coefficients. 

The first calculation of transport coefficients in relativistic scalar theories can be found 



m 



12|. Their calculation is based on Kubo relations, i.e. relations expressing transport 



coefficients in terms of retarded correlation functions between various conserved currents in 
the low momentum, low frequency limit. The correlation functions are directly evaluated 
using finite temperature quantum field theory and thus provide a microscopic calculation of 
transport coefficients; but as is explained in {13! . l^ . their calculation is incomplete. Due 
to the use of resummed propagators (to regularize so-called pinch singularities), an infinite 
number of ladder diagrams must be resummed to get the leading order result. This program 
has been carried out explicitly in ljj, [fjj] for shear and bulk viscosities in scalar theories. The 

TJ3H, 



shear viscosity result has since been reproduced using the real-time formalism 



15 



direct ladder summation in Euclidean space 18] and 2PI effective action methods 



191. The 



2 



results of both shear and bulk viscosities have been reproduced using quantum kinetic field 

□ 

theory derived from the closed-time-path 2PI effective action |20J. 

Order of magnitude estimates of transport coefficients in gauge theories based on phe- 
nomenology appeared more than 20 years ago in |21|, b ut a complete leading order calcula- 



tion in hot gauge theories just came out recently [22|, |23j, |24|, |25(. One of the main reason for 
this is the subtlety of the power counting involved, i.e. which scattering processes should 
be included at leading order. In gauge theories, in addition to 2 —>■ 2 scatterings with a soft 
momentum exchange, the presence of collinear singularities makes the 1 — > 2 scatterings as 



important as the 2 — ^ 2 ones 



26 
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28 



29(. Moreover, interference effects between the 



various collinear emissions must also be taken into account at leading order; this is called 



the 



22 



andau-Pomeranchuk-Migdal (LPM) effect 



23 



24 



29|, 



30|, 



31 



32 



33| . The calculations in Refs. 



25| are based on kinetic theory and consistently include the physics of pinch sin- 



gularities, collinear singularities and the LPM effect. Let us mention that these calculations 
are also very involved technically, since the inclusion of pinch singularities implies the use of 
an integral equation and the inclusion of collinear singularities and the LPM effect implies 
the use of another integral equation embedded in the first one. 

The equivalence between the quantum field theory approach and the kinetic theory ap- 
proach in transport coefficients calculations was shown in [k], [l4] in the case of scalar 
theories; but due to the complications mentioned previously, a similar equivalence has been 
lacking in gauge theories. There exist some attempts at computing transport coefficients 
in hot gauge theories from quantum field theory using different approaches, such as direct 



adder summation in Euclidean space 



18 



34 



351 ] . dynamical renormalization group methods 



36j and 2PI effective action methods 37j, but as far as we know, none of these approaches 



go beyond leading log order accuracy (i.e. with corrections suppressed by 0(g ln^ -1 ))) or 
the large Nf approximation. 

The goal of this paper is to do a leading order calculation of transport coefficients in hot 
gauge theories using purely diagrammatic methods. More precisely, we show the equivalence 
between quantum field theory and kinetic theory for transport coefficient calculations in hot 
gauge theories, thus justifying the effective kinetic theory results of Arnold, Moore and 
Yaffe 23j, |24( from first principles. This nontrivial check is reason enough for doing this 



calculation, but there are other reasons as well that go beyond kinetic theory. For example, 
in cases where the configuration of the field itself is important (e.g. instanton) or unphysical 



3 



particles with indefinite metric appear, then kinetic theory is inappropriate and one must 
rely on quantum field theory. Also, quantum field theory might be the only way of computing 
transport coefficients beyond leading order; according to 23J, doing the calculation using 
kinetic theory would require a whole new machinery. It is not clear if calculations going 
beyond leading order can be converted to a linearized Boltzmann equation. Finally, showing 
that the calculations are grounded and feasible in quantum field theory is important in 
itself. It could give the necessary impetus for other quantum field theory methods (e.g. 
closed-time-path 2PI effective action [20] or the dynamical renormalization group [36]) to 
complete a leading order calculation and provide other insights into the problem. In view 
of the applications mentioned previously, we think it is interesting and important to pursue 
this work. 

In this paper, we restrict ourselves to Quantum Electrodynamics (QED) and electrical 
conductivity; we address the case of shear viscosity in a future paper [38]. The rest of 
the paper is organized as follows. Section [TT1 presents our notation and some background 
material on transport coefficients, both in scalar and gauge theories. Following the work of 
Ref. 



341 ]. Sect. II III presents the derivation of the Ward identity in the limit appropriate for 
transport coefficient calculations. It also presents the constraint on ladder kernels that can 
be obtained from the Ward identity. Power counting arguments are shown in Sect. IIV| in 
order to determine which rungs should be kept in the resummation. The final expression 
for electrical conductivity, including collinear physics and the LPM effect, is presented in 
Sect. |Vl We finally conclude in Sect. ED Technical details of some aspects of the calculations 
are relegated to the appendices. 



II. BACKGROUND MATERIAL 
A. Notation and Useful Formulas 



We present here our notations and some useful formulas that are routinely used through- 
out the analysis. Latin indices run from 1 to 3 and represent space components while Greek 
indices run from to 3 and represent spacetime components. Boldface, normal and capital 
letters denote 3-momenta, 4-momenta and Euclidean 4-momenta, respectively. We use the 
metric convention r)^ = (1,-1,-1,-1). Sums over Matsubara frequencies are written as 



4 



/ (fifr = T S^ P / where P = (ii/ p ,p) and v p = 2mrT (bosons) or u p = (2n + l)vrT 

(fermions) with n an integer. The subscripts B, F attached to a quantity refer to its 
bosonic or fermionic nature (except for self-energies and widths, where we use a special no- 
tation). The subscripts R, I means real or imaginary part and the superscripts ret, adv, cor 
means retarded, advanced or autocorrelation (i.e. average value of the anti-commutator). 
A bar over a quantity means that the gamma matrix structure is explicitly taken out (e.g. 
G"(Jfe) =7"G(Jfe)). 

We give a list of various finite temperature field theory quantities that are used in the 
rest of the paper. We give their explicit expressions in momentum space and not their basic 



definitions in terms of fields (see for example 39, 
our purposes. Free spectral densities are given by 



because the latter are not useful for 



4l|: 



p B (k) = sgn(fc°)27r5((fc i 



0\2 



El) 



p F (k) = 2ir5(k fl - E k )h + (k) +2n5(k° + E k )h-(k) 



(1) 
(2) 



where E k = |k|, h±(k) = (7 =F 7 • k)/2 and k = k/|k|. Note that since we consider systems 
where the temperature is much larger than any other scale, we put m = in the above and 
all subsequent expressions when the momentum of the excitation is hard. From CPT, it can 
be shown that the spectral densities satisfy /?#(— /c°) = — and p F (—k) = pp(k) (in 
the massless limit). At finite temperature, any excitation propagating in a medium has a 
finite lifetime due to numerous collisions with on-shell thermal excitations. The effect o:: 
finite lifetime is to turn the delta functions in Eqs. ((H)-(j2J) into Lorentzians, giving 

1 



this 



pB(k) 
p F {k) 



2E k 



Ik 



Ik 



_(k° 



— E k 



)2 + ( 7fc /2)a (kP + E k y + ( 7fc /2)a 

T k 



;h + (k) 



;h-(k) 



(3) 
(4) 



L(^° - E k f + {T k /2f ^ ' (HP + E k f + (T k /2)* 

The widths are given by 7fe = ^(Jfe = E k )/E k and T k = tr [^} e \k° = E k )\ /2E k , where 
H(k) and S(fc) are the bosonic and fermionic self-energies, respectively. Note that when the 
momentum k is soft, perturbation theory must be re-organized and partial resummation 
of spectral densities is necessary (also called Hard Thermal Loop (HTLs) resummations 



42 



43 



441]). These resummations give rise to screening thermal masses and may also produce 



Landau damping. In gauge theories, HTLs are also essential to obtain gauge invariant results 
(see Sect, lira. 



5 



The time-ordered (or "uncut") propagators can be expressed in terms of the spectral 



densities 



4l|: 



° B ' F{k) = 1 J (2?r) PB/f(UJ) [ kP-u + ie ± ifc°-a;-» e J (5) 
where n B / F (k°) ar e the usual Bose-Einstein or Fermi-Dirac distribution functions. The 



anti time-ordered propagators are just the comp 



Wightman (or "cut" ) propagators are given by |4jJ : 



ex conjugate of the time-ordered ones. 



A% /F (k) = (l±n B/F (k°))p B/F (k) (6) 
A B/F (k) = ±n B/F (k°)p B/F (k) (7) 

The propagators ([5])-(j7j) are the four propagators of the closed-time-path or "1-2" formalism 
46|, with the correspondence G — G n , G* — G 22 , A+ = G 12 and A" = G 21 . Switching 



45 



;o the Keldysh (or r,a) basis, we can also write down the physical functions (see for example 



■/~ira 
lLr B /p 


= iGf /F (k) 


iG a B / F 


= iG°$] F {k) 


■rtrr 


= iG^ /F {k) 



G B/F (k) - A B/F (k) (8) 
G B/F (k) - A+ /F (k) (9) 
A+ /F (k) + A B/F (k) (10) 

The G°^i F is identically zero in the Keldysh basis. Note also that any vertex in this 
basis must involve an odd number of a's (see for example }47|). One can see that 
from the expression of the generating functional in the closed-time-path formalism, Z = 
J V[(p}exp(i J dt J d 3 x (C + J c (f>))- Due to the integration over the closed-time-path, the 
Lagrangian is effectively C = £(<fii) — C(4>2), where 4>i ari d 02 are fields living on the 
time-ordered and anti time-ordered contours respectively. From this we conclude that any 
interaction term is odd under the interchange of <pi and <p2- Switching to the Keldysh basis 
using r = (0x + 02 ) /2 and a = 0i — 2 , we see that any interaction must have an odd 
number of a's, since only a is odd under 0i <-» 2 . As a final remark, to get the explicit 
forms for free or resummed propagators, substitute in the appropriate spectral density (c.f. 
Eqs. ([I])-®). 



The cutting rules we use are the ones that can be found in 
to the zero-temperature ones. The rules are: 



41 



48 



491 ] and are analogous 
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1. Draw all the cut diagrams relevant to the problem considered, where cuts separate the 
unshaded (i.e. "1") and the shaded (i.e. 



"2") regions. 



2. Use the usual Feynman rules for the unshaded region assigning Gb/f{^) to the uncut 
lines. For the shaded region, use the conjugate Feynman rules assigning G* B ^ F (k) to 
the uncut lines. 

3. If the momentum of a cut line crosses from the unshaded to the shaded region, assign 
Ag, F (k). If the momentum of a cut line crosses from the shaded to the unshaded 
region, assign A B , F (k). 

4. Divide by the appropriate symmetry factor and multiply by an overall factor of —i. 
The various shadings given by the cutting rules are not all independent. First, the "vanishing 



of all circlings" relation 
is zero: 



41 



49| says that the sum of all possible cuts of a given diagram 



0,2— Or, 

B/F 







(11) 



a,=l,2 



Equation pip can be rearranged so as to have the same form as a (generalized) optical 
theorem and is thus related to unitarity. The vanishing of all circlings relation is only based 
on the observation that all propagators can be decomposed into a positive frequency part 
and a negative frequency part. Second, the Kubo-Martin-Schwinger (KMS) relations express 
the proportionality between various pairs of shadings. For example, for 2-point functions 
we have [50 ] : 

Gb/f{^) 

For 4-point functions, the relations are 



±e^G% F {k) 



(12) 



51 



52j: 



U B/F 


k 2 , 


k 3 . 


h) 


= ±e-^G 1 B f F (k 1 ,k 2 ,k 3 ,k 4 ) 


U B/F 


k 2 . 


k 3 . 


h) 


= ±e-^G 2 B f F (k u k 2 ,k 3 ,h) 


^B/F 


k 2 . 




h) 


= ±e~^G 2 B 2 } 2 F (k 1 ,k 2 ,k 3 ,k 4 ) 


U B/F 


k 2 . 


k 3 . 


h) 


= ±e-^G 2 iJ F (k h k 2 ,k 3 ,k 4 ) 


r* 2211 (h 

U B/F 


k 2 . 


h, 


h) 


= e-^^G B f F (k u k 2 ,k 3 ,h) 


r* 2121 (h 

U B/F 


k 2 . 


k 3 . 


h) 


= e-^+^G^KhMMM) 


U B/F 


k 2 . 


k 3 . 


h) 


= e-^^G^KhMMM) 



(13) 



where (3 is the inverse temperature and energy-momentum requires k\ + k% + k% + k^ = 0. 
The relations f|T2l) - f|T3l) are only based on the (anti) periodicity in imaginary times of the 
Green's functions and are thus only valid in equilibrium. 



B. Transport Coefficients in Relativistic Scalar Field Theory 

A natural starting point for evaluating transport coefficients using quantum field theory is 
a Kubo type relation. These relations express transport coefficients in terms of long distance 
correlations between conserved currents. Restricting ourselves to electrical conductivity (cr) , 

5J: 



we have [12|, [13j. 



14 



22 



a 



— lim 

6 k°^0, k=0 



d 4 x e*** (j&,x)f(0)) 



<-<! 



(14) 



where j^(x) is the electric current density. Other Kubo relations exist for shear and bulk 



viscosity 



12 



13 



14 



18 



22). Equivalently, one can express these Kubo relations in terms 



of derivatives of the spectral density, itself equal to twice the imaginary part of the corre- 
sponding retarded correlator. Note that the averages are done with respect to an equilibrium 
density matrix, even though transport coefficients are non-equilibrium quantities. 

To compute transport coefficients using these relations, one needs to expand the retarded 
2-point function perturbatively. As explained in details in Refs. [r]. an infinite number 
of ladder diagrams contribute to the transport coefficient, even at lowest order. The reason 
for this is the low frequency limit in the Kubo relation. This limit gives rise to products of 
propagators G(p) with the same momentum (see Fig. [1]). Since finite temperature propa- 
gators possess four poles (one in each quadrant), one faces situations when the integration 
contour is "pinched" between two poles on opposite sides of the real axis in the complex p° 
plane . In equations, we have (E p is the on-shell quasi-particle energy): 



d°p 
W) 



G B {p)G B (p) 



d°p 

W) 
d°p 

(27) \(p° + ie)-E p J V(P' 



Gf(p)GT(p) 
1 



1 



E r 



(15) 



which diverges when e goes to zero. This divergence is symptomatic of the infinite lifetime of 
the excitation. In a medium, excitations constantly suffer collisions from on-shell excitations 
coming from the thermal bath, resulting in a finite lifetime. Thus, the use of resummed 
propagators regularizes these "pinch" singularities, effectively replacing 1/e with l/Tl/(p). 
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p+q 





+ (■■•) 



FIG. 1: Examples of ladder diagrams, where the grey squares represent 4-point functions called 
"rungs". When the external momentum q goes to zero, the two "side rail" propagators have the 
same momentum and produce a "pinch" singularity. 



I F D 
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I 



MFD 



where 



[> 



+ 




FIG. 2: Schematic representation of Eqs. (fT6|) - (fT7|) . The symbols refer directly to the equations: 
2 is an external current insertion, Ai is a rung, T represents a pair of side rail propagators and T> 
is an effective vertex. 



The explicit coupling constants that now appear in the denominator change the power 
counting dramatically and make the resummation of an infinite number of ladder diagrams 

nn 

necessary [13|, H4| • 

This resummation is done by re- writing the infinite sum of ladder diagrams in terms of an 
effective vertex, itself solution to an integral equation. Schematically, we have for electrical 
conductivity: 



cr = — 



6 J (2vr) 4 



T F {k)F(k)V F (k) 



(16) 



V F (k) = l F (k) + / 7^4 ^ } p)V F (p) (17) 

where /C = M.T . See Fig. [2] for a graphical representation of Eqs. f|T6l) - f|T7j) . The symbol 
T represents a pair of "side rail" propagators (note that the "ladder" diagram in Fig. [2] 
is on its side, meaning that the "side rails" are on the top and bottom of the diagram). 
In the limit q — > 0, the two side rail propagators that hook up M. to Dp have the same 
momentum and "pinch", producing a 1/Hi(p) factor. The "rungs" M. are 4-point functions 
that must be of the same order as T~ l but otherwise arbitrary. The effective vertex T> 



9 



encodes the information about the infinite resummation of ladder diagrams. If we close it 
with an external current insertion X (with an additional factor of JF to connect them), we 
get the Kubo relation (1161) . Note that because of the T factor, the transport coefficient 
gets a non-analytic behavior in the coupling constant. This is expected, since transport 
coefficients are roughly proportional to the mean free path and thus inversely proportional 
to the scattering cross section of the processes responsible for transport. 

To write down the appropriate integral equation (TPTj) and solve it are the main tasks of 
transport coefficient calculations. Altough it can become very involved numerically (espe- 
cially for gauge theories), there is no conceptual problem with solving the integral equation. 
On the other hand, to know which rungs contribute at leading order is difficult and requires 

n 

detailed power counting arguments. This is the approach adopted in [13J, where they find 
that a certain (finite) set of rungs are necessary and sufficient to compute the shear and 
bulk viscosities at leading order in a g<p 3 + A0 4 theory. 



C. Additional Complications in Gauge Theories 

Gauge theories are considered to be the most successful theories on the market to describe 
the fundamental interactions of Nature. So for "real" applications (e.g. viscous hydrody- 
namic evolution of the QGP), the need to extend transport coefficient calculations to gauge 
theories is obvious. From experience, we know that gauge theories are more complicated 
than scalar theories, and the present calculations are no exceptions. We make a list of the 
main complications that arise when computing transport coefficients in hot gauge theories. 

Preserving gauge invariance is an important problem when dealing with gauge theories. 
In practice, Ward identities tell us how non-gauge invariant quantities (such as propagators 
and vertices) must be related to each other so as to preserve gauge invariance. We show in 
Sect. |TTl]how Ward identities are used to this effect. 

Another complication, not specific to gauge theories, is related to the presence of different 
species of particles (electrons, photons, etc) and particles with different statistics. The fact 
that there are fermions in the theory means that the tools developed in [3] must be slightly 



modified to take into account the fermionic nature of the particles. These comp 
not major and rather technical in nature. We show in Sect. IVlhow the tools in 
modified in the presence of these complications. 



ications are 



13( must be 
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The distinction between hard (0(T)) and soft (O(eT)) momenta is also particularly 
important in hot gauge theories (e is the electromagnetic coupling constant). In particular, 



it is shown in 



42. 



43 



441 ] that the theory must be partially resummed when soft momenta are 



present to get gauge invariant results. These Hard Thermal Loop (HTL) resummations give 
rise to screening thermal masses and modify the form of the propagators (vertices are also 
modified but this not necessary for our purposes). For numerical purposes (as in 23, 24^). 
the expressions of HTL-resummed propagators are required. In the prese nt p aper, our goal 
is to reproduce the Boltzmann equation of Arnold, Moore and Yaffe 23|, |24J and we never 



use the explicit form of the propagators; the calculations of Sect. [V] are based on very general 
properties that do not depend on HTL resummation, such as even/odd properties of spectral 
densities, KMS conditions and unitarity. In this sense, this technical complication does not 
concern us. 

neories, 



But surely the most important complication comes from the fact that, in gauge t 
transport coefficients are sensitive to soft and collinear physics even at leading order ; 
In scalar theories, explicit^power counting shows that soft momenta and collinear singularities 



23 



24]. 



do not have any effect 13| . As with pinch singularities, collinear singularities come from the 
multiplication of two propagators with momenta p and p + q that are nearly collinear (i.e. 
p-q~0(e 2 T 2 )): 



d°p 
W) 



G B {p)G B {p + q) 



Gf(p)Gt(p + q) 



d°p 

W) 
d°p 

(2^0 \(p° + ie)-E r 
1 

q° + {E p+q - E p ) - 2ie 



1 



(pO + q o _ ie ) _ Ep+q 



(18) 



In the limit q — > (pinch singularities), the expression diverges as 1/e and we get back 
Eq. (fl5|) . In the case where q is nearly on-shell but nonzero and the angle between the 
quasi-particles is small (i.e. 0. 



pq 



~ e), we have E p+q « E p ± |q| and the integral diverges as 
l/6 pq (or 1/e in the perfectly collinear and massless case). These collinear singularities must 
be regulated by including a finite width in the propagators. As with the pinch singularity 
case, the introduction of coupling constants in the denominator changes the power counting 
dramatically. This collinear enhancement affects an infinite class of diagrams: this is the 



well-known LPM effect (see for example 29|, |32j, |33| for a discussion of this effect in the 
context of photon production). Because there are now two types of singularities, there are 
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FIG. 3: Momentum convention for the Ward identity (Eq. (|19|) ). The dotted line stands for an 
insertion of an electric current and the blob represents the amputated fermionic effective vertex 
V F . 

two resummations to do. While in scalar theories there is an infinite number of ladder 
diagrams and a finite number of different rung types, in gauge theories there are an infinite 
number of ladder diagrams (due to pinch singularities) and an infinite number of different 
rung types (due to collinear singularities). We come back to these issues and the subtleties 
of power counting in gauge theories in Sect. IIVI 



III. WARD IDENTITY CONSTRAINTS ON THE LADDER KERNELS 
A. Ward Identities for Charge 

As mentioned in Sect. Ill CI the task of finding which rungs contribute at leading order is 
not easy. The purpose of the present section and the next one is to show that there exists 
a guiding principle that can help us find the structure of the rungs. 

A necessary piece of information in our demonstration is the Ward identity for the effective 
vertex, as derived for the first time for the case of an electric current insertion in Ref. 34]. 



We show here the essential steps. The starting point is the usual Ward identity for charge 
conservation in Euclidean space 

Q^V F (K,K-Q) = G F \K)-G F \K-Q) (19) 

where T>p is a fermionic amputated effective vertex. The momentum convention for Eq. (TT9~j) 
is shown in Fig. [3j To go from Euclidean space to Minkowski space, we need to analytically 
continue K and Q towards real energies. The proper choice here is dictated by the physics of 
pinch singularities and the Kubo relation. To have pinch singularities, one needs the multi- 
plication of two propagators with different boundary conditions, namely G^ t (/c)G^ dv (/c — q) 

12 



or G F iv (k)G T F t (k — q); the Kubo relation requires the evaluation of a retarded current- 
current 2-point function. These two requirements uniquely fix the analytic continuation to 
K — > k° + ie, K — Q — > k° + q° — ie and Q — > q° + 2ie. With this analytic continuation, 
Eq. f[T9"j) becomes 

q,V F (k + ie, k - q - ie) = G F \Jk) - G F \ dv (k - q) (20) 

Taking the necessary limits q° — > and q — > (c.f. Eq. ( fl4|) ) and using C F 1 rct y adv (p) = 
[7°(p° ± iVp/2) — 7 ■ p] (valid when p is nearly on-shell), we get 

lim q.V F (k + ie,k-q-ie) = i-f°T k = 2iZ] et (k) (21) 

The last equality is valid near k° ~ k. This last equation relates the effective vertex of 
the integral equation to the imaginary part of the on-shell retarded self-energy in the limit 
relevant to transport in the case of an electric current insertion. Note that we used resummed 
retarted/advanced propagators because of the need to regularize pinch singularities. As a 



final remark, note that it is possible to repeat the above derivation for a T^ v insertion 38]. 
This is particularly important when calculating shear viscosity from Kubo relations. 

B. Derivation of the Constraint 

In this section, we show that there exists a relation between the kernel of the integral 
equation (TITI) and the imaginary part of the on-shell self-energy. The ladder kernel being a 
4-point function and the self-energy being a 2-point function, the relation between the two 
must necessarily involve functional derivatives with respect to propagators (i.e. "opening" 



of lines in a 
methods 
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eynman diagram). Such relations have been obtained using 2PI effective action 



371 ] . Here we derive similar constraints but starting from a more physical point 
of view, namely the Ward identities of the previous section. 

In the following, we restrict ourselves to electrical conductivity (i.e. only one integral 
equation is involved), but the discussion can be generalized to viscosity with some effort 
(although the result is not as clean cut as for the electrical conductivity) [38j. Starting from 



the integral equation ffTTj) in Euclidean space, 

V F (K, K-Q) = 2%{K, K-Q) + J S^JC(K, P, Q)V F (P, P - Q) (22) 
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our goal is to isolate the kernel and express it in terms of known quantities. Multiplying 
both sides by Q M and using Eq. (fl9|) . we get 

Q,V F {K,K-Q) = Q fl l F (K,K-Q) + J JULm(K,P,Q)[G f (P - Q) - G F {P){2?>) 

where we have separated the kernel into a "rung" part and a "side rail" part, JC(K, P, Q) = 
Ai(K, P,Q)G F (K)G F (K — P). Note that this is a schematic notation and it does not 
reflect correctly the Dirac matrix structure of the integral equation; one must be aware 
of this fact when doing explicit calculations. To make further progress, it is necessary 
to analytically continue towards real energies. This step is delicate, since the sum over 
Matsubara frequencies must be done first. As explained in more detail in Appendix the 
Matsubara sum can be done using the summation formula found in [3] , with M. replaced by 
the spectral representation of a general 4-point function. The physics of the Kubo relation 
and pinch singularities can then be implemented by doing the same analytic continuation 
as for the Ward identities (see Sect. IIII Al) . The result is that the integral equation keeps its 
form in Minkowski space and is similar to the results of Jeon [3] (suitably generalized to 
fermions). Taking the q — > limit and using the Ward identity (J2~T1) . we get 

2iEf(*) = / -0^M(k,p)[Gr( P )~Gf(p)] (24) 

with the understanding that k and p are on-shell in the pinching pole limit. Using the 
definition of the spectral density p F (p) = i(G F t (p) — Gf? v (p)), we can rewrite the integral 
equation as 

2Ef {k) = j J^ M ^ p)p+ ^ p) (25) 

Here S} et (fc) is the full (imaginary) self-energy and Ai is the coherent addition of all possible 
rungs. Since T,} et and pp are known quantities, Eq. (|25|) expresses a constraint on the rung 
part of the ladder kernel M.. To get a more useful expression, we need to invert Eq. (l25l) . 
Let us functionally differentiate both sides with respect to p F (q) 

$mm _ r j^p_ ( sM(k, P ) \ 

S P $(q) ~ J (2-) 4 l 6pt(q) )P^ +M ^ W 
All cut/uncut propagators can be expressed in terms of the spectral density (c.f. Eqs. (jSJ)- 
(0)), thus differentiating with respect to Pp(q) can be interpreted as the opening of a fermion 
line in a Feynman diagram. Such a diagrammatic interpretation is of course not perfect, 

14 



since p~\? is only proportional to fermionic propagators. Some distribution functions are left 
behind when interpreting Eq. fl26|) diagrammatically, but since we are only interested in the 
diagrammatic structure of At, a proportionality relation is sufficient for our purposes. 

Equation f[2T)j) can be further reduced. The key observation is that At is an addition 
of Feynman diagrams made of an arbitrary number of fermion propagators. The following 
relation 

iwf^v^-^ <27) 

is valid for any Feynman diagram A (in particular for Ai and £/ et ), where a is the number 
of fermion propagators in A. In words, it means that the "inverse operator" of a functional 
derivative with respect to a fermion propagator is an integration and a multiplication by a 
fermion propagator (with the same momentum present in the function we differentiate with). 
It is easy to see that Eq. (1271) is true in simple examples (both for bosons and fermions). 

Equation (1271) cannot be applied directly to Eq. (I26I) . because the momenta in the func- 
tional derivative and in the spectral density are not the same. In order to apply the inverse 
operator (!2~Tj) and thus use the special condition that Ai is a Feynman diagram, we need to 
manipulate Eq. (|26|) . Let us functionally differentiate Eq. (l2"oT) a second time, giving 

5M(k,q) 5M(k,r) f d'p ( 5 2 M(k,p) \ 5 2 (2Ef(k)) 

+ = " / 77^4 TZ+TZ^+Ta Pf(p) + T-+7 x,„+/ x (28) 



8p$(r) 5p+(q) ^ J W\5p${r)5p+(q)J , ' Fy " 5p+(r)5p+(g) 

Note that this last equation expresses the "symmetrized" first order functional derivative of 
the rung kernel in terms of second functional derivatives of the rung kernel and self-energy. 
Plugging Eq. (J2SJ) back into Eq. f[2"6T) . we get 

8(2Zf(k)) _ f d*p (5M{k,q)\ p+M+ r dtp { 5 2 (2Zf(k)) \^ p) 



Sp+(q) J (2tt) 4 V Sp+(p) J (2tt)4 \6p+(p)6p$(q). 

f d* P f dH ( 5 2 M(k,l) \ +/ , lAll N 

-J i^AMdm PFm{p)+M[k ' q) (29) 

Written in this form, Eq. (|29|) can be simplified using the inverse operator. For convenience, 
we introduce the following notation 

oo 

M = J2 M « 

a=l 
oo 

E ret = E r /* (30) 

(3=1 
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where the sums over a and (3 are over the number of fermion propagators (note that the 
.M a 's or can be sums of diagrams in themselves). These definitions mean that, in 
general, M. and H] et are "blobs" containing diagrams of all orders; these diagrams can be 
re-organized in terms of their number of fermion propagators. Applying Eqs. (1271) and f[3"Uj) 
to Eq. 0291) and doing some algebra, we obtain 

( w) *» + ^ (If) - «>] " 

(31) 

This last equation is valid component by component, i.e. for each value of the number of 
fermion propagators a. This is a constraint equation, implementing the condition that M. 
is a Feynman diagram and the use of the inverse operator. Plugging Eq. (I3T]) back into the 
first integral equation fl26l) and using Eq. (1301) . we get 

^■>-OT) OT) 

valid for all values of a. The second equality is obtained by summing over a on both sides 
and using Eqs. [30j This is the desired constraint on the ladder kernel. Note that it is 
possible to obtain relations similar to Eq. ( 1321) for shear viscosity [38|]. So even if we restrict 
ourselves to the case of electrical conductivity in QED in this paper, our method is general 
and can in principle be applied to other transport coefficients and other theories. 

C. Interpretation of the Constraint 

As pointed out in Sect. Ill BHH C| one of the hardest tasks in computing transport coeffi- 
cients by diagrammatic methods is to find the right rung kernel in Eq. (fTTj) . In other words, 
the ultimate goal is to find all the appropriate rungs that must be included in Ai so as to 
obtain the transport coefficient at the desired level of accuracy. The usual procedure is to 
consider all possible rungs cut in all possible ways and use power counting arguments to 
show which ones contribute at, say, leading order. This way to proceed is already involved 
in scalar theories; in gauge theories, it would quickly become hard to manage. 

Equation fl32|) provides a way to restrict the number of possible cut rungs that must be 
included in Ai. It says that the full rung kernel is proportional to the functional derivative 
of the imaginary part of the full retarded on-shell self-energy. This result means that any 
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FIG. 4: Rungs obtained by functionnally differentiating with respect to a fermionic spectral density 
the self-energies in the left column (c.f. Eq. [32]) . Power counting arguments show that all these 
rungs contribute at leading order except diagrams (e) and (f), because massless three-body on-shell 
decays are suppressed. All diagrams can be converted into 2^2 scatterings with a soft exchange 
except (h) and (i), which are part of the 1 — * 2 collinear scatterings discussed in 
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20]. 



rung kernel (a 4-point function) is obtained by opening a fermionic line in an imaginary 
self-energy (a 2-point function). The recipe to find Ai is thus to select a certain set of 
imaginary self-energies (using for example an expansion in coupling constants) and open 
the fermion lines in these self-energies to get the rungs. An example of such a procedure 
applied to one-loop self-energies in QED is shown in Fig. |H Since we are dealing with a 
quantum mechanical system, all such openings must be summed coherently. This is a strong 
constraint on the rung kernel. 

A few remarks are in order here. The above constraint implements the physics of transport 
coefficients (i.e. pinch singularities due to the low frequency, low momentum limit) and is 
useful for restricting the number of possible cut rungs: it is a necessary condition to get 
leading order rungs. In this case, power counting is used to verify if the rungs obtained 
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from the constraint contribute at leading order. For electrical conductivity, this verification 
is done in Sect. IIV A[ If other singularities are present, then the constraint is not sufficient 
and must be supplemented with additional power counting arguments. This is the case in 
gauge theories, where collinear singularities make rungs with number changing processes as 
important as those without these processes. These additional power counting arguments are 
presented in Sect. IIVBI 

The constraint equation fl32|) expresses rungs in terms of functional derivatives of self- 
energies. Since a few self-energies can give rise to many rungs, it could be more convenient to 
work directly on the self-energies in some cases. This is particularly relevant when, because 
of the presence of collinear singularities, an infinite number of rungs are needed at leading 
order. Such an analysis could rapidly become unmanageable; as shown in Sect. IIVB[ using 
power counting to carefully select the relevant self-energies before opening them makes the 
analysis easier. 

The constraint ( l32|) is also essential to preserve gauge invariance in our formalism. To 
see that more clearly, let us compare our formalism to the 2PI formalism. Expressions 
similar to Eq. ( 132]) were obtained using 2PI effective actions methods Q , with notable 
differences. In 2PI methods, the "constraint" is in coordinate space and comes naturally 
from standard functional relations. The important point is that the kernel M.2PI of the 
integral equation is given by the functional derivative of the self-energy with respect to 
resummed propagators, where the self-energy is itself given by the functional derivative of 
all amputated 2PI diagrams ^2) with respect to resummed propagators: 

To obtain Eq. fl5S|) . no reference is made to the low frequency, low momentum limit or pinch 



singularities. In that sense, the "constraint" in 
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37j is too general for the problem at 



land. We also mention that there seem to be gauge invariance issues with 2PI methods 



54 



55, 



56j|. Briefly, this is because the 2PI effective action methods involve the use of 



resummed propagators, i.e. only a certain class of topologies are resummed. On the other 
hand, gauge theories require Ward identities to be satisfied. Since Ward identities require 
cancellations between different classes of topologies (resummed propagators and vertices), 
they are necessarily violated in the 2PI formalism. Note that it is shown in [57| that a 
complete, self-consistent system of equations for the dynamics of QED or QCD can be 
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obained using higher order effective actions, but it has never been implemented in practice 
for transport coefficients. 

In contrast, our method does not make any reference to 2PI effective actions and starts 
directly from symmetry principles, the Ward identity being the expression of these symme- 
tries for quantum mechanical amplitudes. The physical limit that is used to obtain transport 
coefficients (q — > 0) is explicitly implemented in our constraint (Eq. (|32|) ). In addition, the 
constraint contains, by construction, exactly those diagrams required to produce a gauge 
invariant result for the electrical conductivity in QED; in other words, it tells us which 
self-energy must be resummed in the side rail propagators for each rung present in the ker- 
nel so as to keep everything transverse. This constraint is not as general as the constraint 
obtained from 2PI methods (compare Eqs. f[52"j) and (1551) ) or from general considerations of 
gauge invariance, but it is more powerful in the sense that it contains the relevant physics 
to calculate the electrical conductivity in QED. 

IV. POWER COUNTING 

Equation (132]) gives us an infinite number of self-energies from which rungs can be ob- 
tained by "opening" lines. To make progress, we need a selection criterion based on power 
counting to isolate the rungs that contribute at leading order. Naively, one could do an ex- 
pansion in coupling constants and keep only the lowest order diagrams, but this procedure 
turns out not to be sufficient. As pointed out in Sect. Ill C] due to the presence of collinear 
singularities, there is a restricted but infinite class of diagrams that must be resummed to 
get leading order results. The goal of this section is to present power counting arguments to 
obtain the leading order electric conductivity. For the purpose of this section, we divide the 
rungs into two categories: the ones containing collinear singularities and the ones that do 
not. In kinetic theory language, these correspond tol + iV— > 2 + iV and 2 — > 2 scatterings, 
respectively (we come back to this last point in Sect. M) . 

A. Power Counting Without Collinear Singularity 

We first consider the case of no collinear singularity. According to Eq. ([52]) . the rungs 
are obtained by opening fermion lines in the imaginary part of the retared self-energy of the 
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electron. At one loop, the imaginary part is zero since an on-shell massless excitation cannot 
decay into two on-shell massless excitations. It is thus necessary to go to two loops for a 
leading order result. The two-loop imaginary retarded self-energies are shown in Fig. HJ with 
their corresponding rungs. Note that there are many more cuts that correspond to imaginary 
two-loop self-energies. We consider these other cuts when writing down the 4x4 matrix 
integral equation (c.f. (T40l ); for the moment, we are only interested in the rung topology. 
Let us check the power counting size of each rung. 

Consider rung (a) in Fig. HI reproduced in Fig. [5] with momentum labels. The expression 
for the rung is: 

M [a) = -*e 4 f A~(k - Z)A+(p - l)G* B (l)G B (l) (34) 

where for clarity we omitted Dirac matrices and Lorentz indices (irrelevant for non-collinear 
power counting). If all momenta are hard, then the rung is 0(e 4 ) (this is true for all rungs 
considered here). In the case of a soft bosonic exchange / ~ O(eT), the story is different. 
Each bosonic propagator is now 0(e~ 2 T- 2 ) (note that there is no contribution from the 
Bose-Einstein distribution since I is off-shell) and there is a phase space suppression of 
e 2 (i.e. two integrals are killed by the delta functions inside the cut propagators, leaving 
d 2 l ~ 0(e 2 T 2 )). Collecting all powers, we get that rung (a) is 0(e 2 ) for a soft exchange. 
This is the usual Coulomb divergence in scattering theory and is the dominant part of the 
integral. Rungs (b) and (c) also contain a Coulomb divergence; the power counting is done 
in the same way and gives 0(e 2 ). 

The structure of rung (d) is very similar to rung (a), with the role of bosons and fermions 
interchanged. The only difference in the power counting comes from the fact that the 
exchanged soft excitation is now a fermion. Since a (soft) fermion propagator is 0(e _1 T _1 ) 
(instead of 0(e~ 2 T -2 ) for a soft boson propagator), rung (d) is suppressed by two additional 
powers of e compared to rung (a), giving 0(e 4 ). Rungs (e) and (f) are necessarily higher 
order than the others and should not be considered. This is due to the presence of three- 
point vertices with their three legs on-shell and massless, corresponding to kinematically 
forbidden or highly suppressed processes. 

The power counting of rung (g) in Fig. H] is slightly more involved. The expression for 
the rung is (the momentum convention is shown in Fig. [5]): 

/ J0F M V = I I ^ A B( l )& B (k-l-p)G* F (k-l)G F ( P + l) (35) 
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FIG. 5: Momentum labels used to do the power counting of rungs (a) and (g) (c.f. Fig. 2]). 



where we explicitly write the integration over p coming from the integral equation (c.f. 
Eq. (TITj) ). We consider the case where the side rail momenta k and p are both hard and 
nearly on-shell (i.e. k 2 ~ p 2 ~ 0(e 2 T 2 )) while the loop momentum I is soft. In such a 
kinematical regime, the dominant contribution comes from when the propagor momenta 
satisfy {k — l) 2 ~ (p + l) 2 ~ 0(e 2 T 2 ). For these conditions to be satisfied, the angles 6ki and 
9 p i must be 0(e), implying an e 2 suppression in both phase space integrations (this is because 
d 4 p = dp \p\ 2 sm8 d\p\ dOdcf) ~ 0(e 2 ) when 6 ~ 0(e)). Another way to see the phase space 
suppression in dH is to notice that the two delta functions inside the cut propagators kill two 
integrals over /, leaving only d 2 l ~ 0(e 2 ) when I is soft. Accounting for all powers of e, we 
have an e 4 coming from the four explicit vertices, e 2 x e 2 from the phase space integrations 
and e~ 2 x e~ 2 coming from the two fermion propagators (note that the fermion propagators 
are not 0(e^ 1 T^ 1 ) because their upstair momenta are hard). Rung (g) is thus 0(e 4 ). The 
power counting and size of rungs (j), (k) and (1) are similar to rung (g). Finally, rungs (h) 
and (i) are suppressed except in the very special kinematical regime where the electron is 
collinear to the exchanged photon; they are part of the "collinearly singular" rungs and will 
be considered shortly. 

From the above power counting, we get that the rungs with a Coulomb divergence (i.e. 
(a), (b), (c)) are 0(e 2 ) and the rungs without this divergence (i.e. (d), (g), (j), (k), (1)) are 
0(e 4 ). Naively, one would expect only rungs (a)-(c) to contribute at leading order, but it 
turns out that there exists another suppression mechanism that is only effective for these 
rungs. To understand this suppression, let's look at the 2 — > 2 collision term of the linearized 



Boltzmann equation of Arnold, Moore and Yaffe 



221 ] (we show in Sect. IVlhow to obtain this 
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241 ]. The labels u-n represent the species (fermion or photon), flavor and helicity 
of the excitations (here u = electron). A tedious calculation shows the equivalence between the 
scattering processes and the rungs of Fig. SJ The correspondence goes as follows (the letters refer 
to diagrams in Figs. H] and [U respectively): (a) = (b) 2 and (d) 2 , (b) = (e) 2 and (f) 2 , (c) = (a) 2 and 
(c) 2 , (d) = (g) 2 and (h) 2 ,(g) = (gJCh), (j) = (e)(f), (i) = (a)(b) or (c)(d), (j) = (c)(d) or (a)(b). 



equation, c.f. Eq. (H9j) ): 

f,s,h 



C u {k) 



1 



E 



d 3 p 



d 3 l 



'37' 



d 3 l 



(2tt) 3 2E p (2tt) 3 2Ei {2n) 3 2E v 



\M„ 



(k, P] l,l')\ 2 (27rY5\k + p-l-l') 



x n u {k)rf{p){\ ± n m (l))(l ± n n {l')) [ X u (k) + X v (p) - X m (l) - X n (l')] 



(36) 



where the sum is over flavors (/), species (s) and helicities (h), the M uvrn nS are the 2^2 
scatterings relevant for the calculation and the x' s represent small deviations away from 
equilibrium distribution functions. The key observation to understand the suppression is 
that rungs (a), (b) and (c) correspond to 2 — > 2 processes for which both incident excitations 
undergo a soft scattering (k — p = q, where q is O(eT)) without changing their species types 
(u = m and v — n); see Fig.[6]for the correspondence between rungs and scattering processes. 
In such a case, there is a partial cancellation between the x's in the collision term 22]: 
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X u (k)- X u (p) = - q .V X u (k) + 0(g 2 ) ~ 0(eT) + 0(q 2 ) (37) 



Since in the Kubo relation ( |T6i) the term [x w (k) + x"(p) ~ X m (\) ~ X™(1')] is squared, the 
overall suppression is e 2 . Taking into account this additional suppression, rungs (a), (b) and 
(c) are 0(e 4 ). 

In summary, rungs (a)-(c), (d), (g) and (j)-(l) in Fig. [Hall contribute at leading order for 
electric conductivity, although the contribution of rungs (a) and (b) are effectively zero due 
to Furry's theorem. This completes the power counting for leading order rungs corresponding 
to 2 — ► 2 scattering processes in QED. Note that the Ward-like identity constraint is also 
applicable to scalar theories; it can easily be shown that takin g al l the two-loop imaginary 
self-energies of g0 3 + A0 4 gives the necessary leading order rungs 13(. We can conjecture that 
it is necessary and sufficient to use the Ward identity constraint and an expansion in coupling 
constants to get all rungs of a given accuracy that do not depend on another singularity (such 
as collinear singularities). On the other hand, the constraint only implements the physics of 
pinch singularities and thus cannot give precise information about rungs containing collinear 
singularities. 



B. Power Counting With Collinear Singularities 

Let us now come back to the problem of collinear singularities. As explained in Sect. Ill C| 
collinear singularities appear when two propagators G(k) and G(k + p) nearly pinch, i.e. 
when the angle between k and p is Q^ v ~ 0(e). Thus we must look for cut rungs containing 
number changing processes and pairs of nearly pinching propagators. By inspection, we find 
that the only rung topologies satisfying those criteria are the ones with vertex corrections, 
such as rungs (h) and (i) in Fig. HI 

Let us check the power counting size of rung (h), reproduced with momentum labels in 
Fig. [7] (rung (i) is done in a similar way). The expression corresponding to rung (h) is: 

X'fG' F (k + l)'t S &?(k + P + trfVpip)'? (38) 

where we explicitly write the integration over p and the effective vertex coming from the 
integral equation (c.f. Eq. (jTTj) ). The integral over dl° is dominated by the kinematical range 



23 



1 




FIG. 7: Momentum labels used to do the power counting of rung (h) (c.f. Fig. 

where the two fermionic propagators k + I and k + p + I nearly pinch. In such a regime, 
expression ( 1381) becomes: 

d 4 p if d^p f d 3 l 



\p° + (Ek+i — E k+ i +p ) — |(r fc+ z + T k+p+ i) J 

The external legs are hard and nearly on-shell due to pinch singularities, i.e. k ~ (k+p) ~ T 
and fc 2 ~ (fc+p) 2 ~ 0(e 2 T 2 ). We also assume that the exchange photon p is hard (otherwise 
the rung would be subleading due to too much phase space suppression); it is also on- 
shell since it is cut. The nearly pinching conditions are given by (k + Z) 2 ~ 0(e 2 T 2 ) and 
{k + p + l) 2 ~ 0(e 2 T 2 ). This is equivalent to the statement that I 2 , (k-l) and (p-l) must be 
0(e 2 T 2 ) or that the electron on the upper rail must be collinear with the exchange photon. 
In particular, it implies that 9 p i ~ 0(e), thus restricting the phase space of p to an 0(e 2 ) 
region since d 4 p = dp \p\ 2 sin8 d\p\ d9d<p ~ 0(e 2 ) when 9 ~ 0(e). 

Up to now, we have not specified the vertex correction momentum I: it can be either 
hard or soft. Let us first consider the case when I is hard. In this momentum regime, the 
power counting size of rung (h) is e 4 x e 2 x e 2 x [A^JZ)] x py a ($ +/)7 <5 (^ + j4 +/)7 /3 ] x e~ 2 , 
where the e 4 comes from the four explicit vertices, e 2 x e 2 from the (small angle) restriction 
on the phase space of p and /, and e~ 2 from the nearly pinching propagators. Since / is 
hard, the bosonic cut propagator is 0(e°) and is proportional to g a p. In this case, the Dirac 
structure reduces to a scalar product and generates an additional e 2 suppression. Collecting 
all powers of e, we get that rung (h) is 0(e 8 ) when / is hard. 
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The situation is different when I is soft. In this case, the power counting size of rung 
(h) is e 4 x e 2 x e 3 x [Aj a/3 (Z)] x [y a (^ + /)7 5 (^ + i> + /)7 /9 ] x e~ 2 , where the e 4 comes 
from the four explicit vertices, e 2 from the restriction on the phase space of p, e 3 from 
the soft integration over I and e~ 2 from the nearly pinching propagators. Since I is soft, 
the cut propagator is HTL resummed and is no longer proportional to g a p. There is thus 
no additional e 2 suppression coming from the Dirac structure when / is soft. The power 
counting size of A^ HTL (7) ~ (1 + iib(1 ))pb htl(0 is also variable and depends on the 
momentum flowing through it. For soft spacelike momenta, Landau damping gives rise to 
an 0(e 2 ) imaginary self-energy. In this situation, Pshtl(0 ~ ^-f(0/(^ 2 + ^(0) 2 ~ e ~ 2 
and A^ HTL (7) ~ e~ 3 . Collecting all powers, we find that rung (h) is 0(e 4 ) in the collinear 
regime. On the other hand, for null or soft timelike momenta, ~ 0(e 3 ) (no Landau 

damping) and Pb htl(0 ~ e -1 , making the rung subleading. 

More generally, rungs obtained by opening lines on self-energies with any number of vertex 
corrections (as shown in Fig. [8]) could be leading order (see below for an explicit example of 
power counting with two self-energy corrections). The only restriction is that one should not 
open a pinching propagator; since pinch singularities occur for pairs of propagators, opening 
(removing) one of them just cancels the effect and would thus make the rung subleading. 
With this restriction in mind, opening the self-energies in Fig. [8] results in the rungs of 
Fig. [9j Before going further, let us verify the power counting size of the rung in Fig. [7] with 
an additional vertex correction either on the upper or lower rail. It modifies Eq. (|38|) in the 
following way: it adds two explicit vertices (e 2 ), one integral over soft 3-momenta (e 3 ), two 
pinching propagators (e -2 ) and one cut propagator with soft spacelike momentum (e -3 ). 
Collecting together the newly introduced powers of e, we see that adding a vertex correction 
results in an 0(1) correction. This is true for any number of vertex corrections. 

There is one further simplification that can be used to reduce the number of rungs in 
Fig. [91 To see this simplification, rungs must be expressed in the Keldysh basis. In this 
basis, diagrams are expressed in terms of retarded, advanced and correlation functions (c.f. 
Eq. (jSJ)) and it is thus easier to isolate pinching contributions. This imposes constraints 
on the r,a structure of rungs. In particular, the 4-point rungs must be Ai(k, —k, —p,p) = 
■Maarr(k, —k, —p,p) or M.rraa.(k, —k,—p, p) in order for the side rails to pinch (a similar 
argument is used by the authors of [la . Il6l ] to simplify their integral equation) . To have nearly 
pinching contributions from the collinear electron, we also need an alternation of retarded 
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FIG. 8: Self-energies considered in our leading order electrical conductivity calculation. Soft 
photons are indicated by thick lines; all other excitations are hard. In each diagram, the hard 
incoming electron is collinear to the hard photon. Any number of soft vertex corrections must 
be included due to the presence of nearly pinching pairs of progagators; this is the diagrammatic 
depiction of the LPM effect. The propagators indicated by vertical arrows are the only one that 
can be opened without breaking a pair of nearly pinching propagators. 

(ra) and advanced (ar) propagators on each side of the hard collinear photon. Figure [TU] 
shows typical rungs of Fig. [9] in the Keldysh basis submitted to the above constraints. As 
can be seen from the figure, rungs with vertex corrections on both the upper and lower rails 
always contain an aa propagator and are thus zero. Other r,a combinations are possible, 
but they do not pinch and are thus sub leading. In contrast, rungs with vertex corrections 
only on the upper or lower rail can have a pair of pinching propagators and are thus leading 
order. 

This infinite number of rungs with arbitrary number of vertex corrections must be re- 
summed using another integral equation. This is the diagrammatic implementation of the 



LPM effect 
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331 ] . Physically, the LPM effect arises because of coherence effects in 
collinear photon emission/absorption. In the process, the electron can suffer an arbitrary 
number of collisions with soft photons coming from the medium; since the formation time 
of the emitted/absorbed photon is of the same order as the mean free time between soft 
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FIG. 9: Possible rung topologiess containing 1 — ► 2 collinear processes obtained by opening (al- 
lowed) fermion lines in the self-energies of Fig. Soft photons are indicated by thick lines; all 
other excitations are hard. The leading order rungs are in the first row and column; all rungs with 
at least one self-energy correction on both the upper and lower rails are subleading because of the 
impossibility of having upper and lower pairs of propagators pinch at the same time. 

scatterings, then all the processes must be added coherently. We will see more clearly in 
Sect. fVlhow the rungs in Fig. [9] are related to these 1 + iV — > 2 + iV collinear processes. 



V. DERIVATION OF THE INTEGRAL EQUATION 

Sections IIIIBI and JV] show what "rungs" , in the realm of all possible rungs, should be 
resummed in order to obtain the electrical conductivity at leading order. A visual summary 
of the necessary rungs is shown in Fig. [TTJ These rungs can be resummed using the integral 
equation i fPTj) in order to get the electrical conductivity. As explained in Sect. II V[ the rungs 
fall in two categories: those corresponding to 2 — > 2 scatterings (i.e. without collinear 
singularity) and those corresponding to 1 — >• 2 collinear scatterings (i.e. containing collinear 
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ra propagator... ^> one must be ar 

a r a P r r 




FIG. 10: Typical rung topologies for collinear processes expressed in the Keldysh basis. To have 
a leading order contribution, the rungs must be A4 aarr (k, —k, —p,p) or A4 rraa (k, —k, —p,p) and 
must contain an alternation of ra and ar propagators. Taking into account those constraints, it is 
easy to see that for all r,a structures, there is always an aa propagator in topology (b); it is thus 
suppressed compared to topology (a). 

singularities). Note that there is a finite/infinite number of rung types corresponding to 
2 — s- 2/1 — > 2 scatterings, as stated in Sect. Ill CI 



equat 
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ion and show its equivalence 



24} obtained using an effective 



The goal now is to write down the appropriate integra 
to the leading order results of Arnold, Moore and Yaffe 
kinetic theory. The plan is to separate the analysis of collinear and non-collinear rungs in 
equation (TPTj) according to K — (A4 + M)J r ) where Af and M. correspond to the collinear 
and non-collinear rungs, respectively. 



A. Integral Equation Without Collinear Singularity (A4) 

For simplicity, let us first consider the case without any collinear physics. The diagrams 
that need to be resummed are shown in Fig. [Ill To write down the integral equation, we 
follow the general diagrammatic method of Ref. [isj. In this method, all possible ways of 
cutting the effective vertex and the kernel are considered, resulting in a 4 x 4 integral matrix 
equation: 

= + / 70i M(k,p)F( P )T>»(p) (40) 

We use boldface letters here to emphasize the fact that T> and X are 4 component column 
vectors and 7C = is a 4 x 4 matrix. See Fig. [12] for the graphical representation of 
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FIG. 11: Diagrammatic summary of our leading order calculation of electrical conductivity in hot 
QED. The basic integral equation (c.f. Eq. (|17p ) is due to the usual pinch singularities; its solution 
(represented by a grey half circle) must be substituted in the initial Kubo relation to get the 
conductivity. All the rungs included in the kernel of the basic integral equation (represented by a 
grey rectangle) are shown. These rungs can be divided in two categories: those corresponding to 
2 — > 2 scatterings (first and second lines) and those corresponding to 1 — > 2 collinear scatterings 
(third line). The collinear rungs represent an infinite number of vertex corrections (represented 
by a grey triangle) that are resummed using the collinear integral equation (c.f. Eq. (|50p ). Also 
shown are the self-energies that must be resummed in the side rail propagators in order to preserve 
gauge invariance, as dictated by Eq. ([32]) . 
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FIG. 12: Graphical representation of the integral equation (|40p . The inhomogeneous term T(p) is 
not shown explicitly. It is a column vector with only the first component being nonzero, since it is 
an operator insertion and is thus point-like (i.e. it can only be in the shade or out of the shade). 



this equation. We again emphasize that, here and in the following, the order of the various 
fermionic components is not respected; for example, the effective vertex T> F (p) has a Dirac 
structure and should be sandwiched between the two fermionic propagators contained in 
ZF(p) , something that is not apparent from the present notation. Only explicit calculations 
show that the Dirac structure all works out. The matrix Al corresponds to the "rungs" of 
the ladder diagrams and J- corresponds to the "side rails" . The explicit form of J- is: 

/ G F (P)G F ( P ) A F (p)A F (p) G F (p)A F (p) A + F ( P )G F (p) \ 

A F (p)A F (p) G* F (P)G* F (P) & F (p)G* F (p) G F (p)A F (p) 

G F (p)A F (p) A F (p)G F (p) G F ( P )G* F (P) A + F (p)A F (p) 

\ A F (p)G F (p) G F (p)A-(p) A-(p)A F (p) G F (p)G F (p) ) 



Hp) 



i 



-2nAl 



2fl-n)Al 



(1 - 2n)A F 

~2nA F -2nA F (1 - 2n)A F 

(l-2n)A+ (l-2n)A+ ((1 — n)A F — nA F ) 
\ (1 - 2n)A F (1 - 2n)A F —2nA F 



(1 - 2n)A+ 
(1 - 2n)A + F 
2(1 -n)A 
((l-n)A3 



F 

nA~, 



(41) 



/ 



where we have used n = n F (p°) and A^ = A^(p°) and taken the pinch limit (q — > 0) in the 
last line. In matrix form, Eq. (1401) has a complicated structure and is quite cumbersome. In 



the scalar case, it is shown in 



13] that the 4x4 integral matrix equation can be reduced to 



a one- dimensional integral equation using a special transformation. This is expected, since 
one could start directly from the (one-dimensional) Euclidean integral equation, do the 
Matsubara sum and do the analytic continuation, resulting in a one-dimensional integral 
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equation for real energies that should be equivalent to Eq. (140]) (see for example [l8( for 
such an approach). Thus following Ref. [l3|, we decompose the side rail matrix f]4T]) into 
outer products: 



F(p) = w{p)u T {p) + h(p)f(p) (42) 



where the vectors w(p), u(p), h(p) and j(p) are given by: 



wT{p) = 2n F (p°)A F (p) / i i (1 _ e ^) (1 _ e - gp °) 



F V " " 4 4 



J T (p) = (o 1 e^°) (43) 

This decomposition into outer products helps in simplifying Eq. (14*0]) . To see that, substitute 
Eq. (1421) into Eq. (HOI) and iterate a few times. Multiplying the result by u T (k) on both sides 
and noting that j T (p)X(p) = 0, u T {p)h{p) = 0, j T (p)w(p) = 0, j T (p)A4(p, /) oc j T (p) and 
Ai(p, l)h{l) oc 7t(p), it is possible to show that the term h(p)j T (p) in the decomposition (142]) 
gives no contribution when iterating. The integral equation thus becomes: 

u T (k)T> F (k) = u T (k)X F (k) + J -0^[u T (k)M(k,p)w(p)][u T (p)'D F (p)] (44) 

We can relabel the reduced effective vertex by T> F (k) = u T (k)T^ F (k). Since only the 
first component of X F {k) is nonzero, we have u T (k)T F (k) = —X F {k). If we define 
M! = u T (k)A / l(k,p)w(p) as our reduced kernel, the resulting integral equation becomes: 

V F (k) = -l F (k) + J 0^M'(Kp)V F ( P ) (45) 

This last equation can be further reduced. Using the condition of unitarity ( fTT]) to get 
rid of some Ai components, the KMS relations for 4-point functions (I13p and the relations 
M 02 = M{ 2 , M20 = Mix, = M*22 and M 32 = M* Z2 (these last relations must be 

checked explicitly for each rung), a tedious calculation shows that Eq. (145]) reduces to: 

d 4 p 



V l F (k) = -T F (k) + J 



(2tt)« 

24 (p) 



n F (p°)(e^ k ° + 1) [M 22 (k,p) + e^°M 32 (k,p)) (e*° + l)A+(p)J ^^46) 



p 
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where the A4i/s = 0, ...,3) refer to the matrix components of (correspond to the 
different ways of cutting the rung kernel, see Fig. fl~2]) . we defined A^(p) = ^>Ap(p) and we 
used the fact that T> oc 7*. Equation f)46p is the fermionic version of Eq. (1A15I) first obtained 
in [lsj] for scalar theories. The fact that only "completely cut" components of M. appear in 
Eq. (jinD is consistent with the Ward identity constraint (1321) . However, it is still necessary to 
go through all the analysis to get the right proportionality factor, since the constraint (|32|) 
is a proportionality relation when interpreted in terms of diagrams. The reason why it is 
important to consider all possible cuts of M. is because there are many ways to get an 
imaginary part at finite temperature; for example, we have 21m E 11 = E 11 + X 22 = X 12 + X 21 
for any (time-ordered) two-point function (a consequence of unitarity). 

As can be seen from Fig. HJ the rungs are all made of 4-point functions with two external 
vertices in the shade and two out of the shade. Since cut propagators represent nearly 
on-shell thermal quasi-particles, the rungs can be naturally interpreted as 2 — > 2 scattering 
processes. A tedious calculation shows that the sum of all the leading order rungs in Fig. |4] 
can be converted into the square of a scattering matrix, where the scattering processes are 
given in Fig. EJ It is easy to see that diagrammatically, starting from the scattering processes 
(see caption of Fig. [6] for details). In its present form, the effective vertex in Eq. fj46l) is a 
matrix in spinor space; for the Kubo relation (fl6l) to be a scalar equation, we want T> F (k) to 
be a vector in spinor space. Multiplying Eq. fj46|) from the left with u s (k) (or equivalently 
v s (k), see 18] for their definition) and defining D F (k) = u s (k)V F (k), we get: 

d A p d 4 l d A V 



om = -im + j ^^(^( 2 ^ (4) ( fc +^-^0(e^° + i) 



x 



1 f,s,h 

2 £ |M„(fc,p;/,0| 2 A+(Z)A+(0A+(-p) 



D, 771,71 



dHp) 



(47) 



where the M uvm nS are the 2^2 scattering processes shown in Fig. [6j The sum is over 
flavors (f), species (s) and helicities (h). Note that for electrical conductivity, v can only be 
an electron or a positron, otherwise the resulting rung would not connect with the fermionic 



effective vertex (see Fig. EI 



equation obtained in 
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Equation (T47j) is formally equivalent to the linearized Boltzmann 



241 ]. with all leading order 2^2 scatterings included. ' 



"o show 
3 and 



this more explicitly, we define the deviation from equilibrium x'(^) = D l (k)/Tk 
substitute it into Eq. fT4T|) . We also use the delta functions contained in the cut propagators 
to put the integrated momenta on-shell. We remark here that at finite temperature, both 
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positive and negative energies appear in the cut lines. The resulting integral equation is: 

d 3 p dH dH' 



f F (k) = T kX l F {k) + / 



(27T) 4 5 (4) (fc+p-/-/') 



X 



(2tt) 3 (2tt) 3 (2tt) 3 
iV^| M (hnl ,w n v (p°)(l±n m (l°))(l±n n (l' )) 

2 ^ \Muvmn{K, P ,l,l)\ _ 



UUP) - Xln(l) - X»(0) 

(48) 



with the understanding that X7 = (i- e - there is no bosonic effective vertex). Using the 
fact that Ylivmn \M uvmn (k,p;l,l')\ 2 is invariant under the change of labels p <-> — / <-> /' 
and the relation = (1 + e _ ^ fc °)S + (A;) (where S + (fc) is a Wightman self-energy that can 
be expressed in terms of \ Y^vmn \ M uvmn(k,p; I, l')\ 2 ), we finally get: 



n F {k°))P F {k) 



d 3 p dH dH' 



f,s,h 



(2tt) 3 (2tt) 3 (2tt) 3 
xn F (k°)n v (p°)(l ± n m (l°))(l ± n n (l' )) 

x ix\k) + x i(p) - xUi) - xun] 



(27r) 4 5^(k + p-l-l') \M uvmn (k,p;l,l') 



(49) 



This last equation is identical to the one obtained by Arnold, Moore and Yaffe 23|, |24| using 
effective kinetic theory. 



B. Integral Equation With Collinear Singularities (M) 

The case with collinear singularities can be studied with the tools developed in the pre- 
vious section. The only difference comes from the cut rungs that must be included in the 
4x4 integral equation. In Fig. [121 the rung matrix is made of 4-point functions cut in all 
possible ways. For collinear singularities, power counting shows that an infinite number of 
3-point functions, corresponding to an infinite number of vertex corrections (see Fig. [TBI) , is 
needed at leading order. Fortunately, this infinite number of 3-point functions has a simple 
structure and is thus manageable. In analogy to the "usual" transport coefficient calcula- 
tion, we can define an effective vertex V F (k,p) that satisfies the following integral equation: 

j^KoA p, q)F(k, p, q)V F (p, q) (50) 

where 7V C oii is a rung with a single soft photon exchange and the external photon is collinear 
with the electron. See Fig. [T3] for an illustration of this integral equation. This resummation 
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FIG. 13: Leading order rungs necessary for the calculation of electrical conductivity. This infinite 
series of rungs can be resummed into an effective vertex Vp(k,p) using the integral equation (|50p . 






FIG. 14: Illustration of the integral equation (|50p . The kernel Af co n is made of only a single soft 
photon exchange. 

is relevant for photon production including the LPM effect and is done in great detail in 
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581 ] . Instead of using the closed-time-path or "1-2" formalism, they use the r/a formalism 
and are able to put the integral equation in a form convenient for numerical purposes. 

Coming back to our initial problem, we need to find the components of the 4x4 rung 
matrix. In the case of 4-point functions, the 16 entries of the rung matrix are given by the 
different ways of cutting a rung with 4 external vertices. The collinear rungs considered 
in Fig. [13] are 3-point functions with one internal vertex. The cuts coming from the three 
external vertices fill 2 3 = 8 entries of the rung matrix and each entry has two cuts due to 
the internal vertex. Since the rung matrix is not completely filled, the reduction procedure 
used in Sect. IV Al to get Eq. (j4"B"l) could be spoiled and the whole technique would break 
down. Fortunately this is not a problem, because each of the two topologies considered in 
Fig. [13] fill different parts of the rung matrix. It is easy to check that the discussion preceding 
Eq. PS"]) is still valid in the present case. 

Applying the reduction procedure of the preceding section to the topologies shown in 
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a) 



b) 





c) 



d) 





FIG. 15: Cut rungs left after the reduction procedure of Sect. IV Al The cuts in the left column 
are included in A/22; the others were already included in the non-collinear analysis. 



Fig. 



we thus get an equation similar to Eq. fl46|) : 

d A p 



V F (k) = -T F {k) + 



(2tt)< 



V F (p) 



r r 



(51) 



n F (p°)(e^ k ° + l)N 22 {k,p){e^° + l)A+(p) 

where the A/i/s = 0, ...,3) refer to the matrix components of A/ (correspond to the 
different ways of cutting the rung kernel, see Fig. [P2j) . Note that A32 is automatically zero 
for a 3-point function; this is why it is not present in the above equation. The only cuts 
left after the reduction procedure are shown in Fig. [H3 Of the four cuts shown, only the 
ones in the left column are included in A/22; the ones in the right column are part of the 
2 —>■ 2 scatterings considered in Sect. IV Al and are left out of the collinear analysis. The 
collinear rungs (a) and (c) can be written as A/22 (a) = [ — ^VfJAJ (A; — p)[— iVjt^ 2 ]* an d 
•A/22 (c) = [ — iVi^^B au/(^ ~ P)[~ *Vf]*' i- e - as a multiplication of a bare 1 — > 2 vertex V^ 2 
and an effective vertex V F . Note that rungs (a) and (c) are complex conjugate of each other; 
this is a consequence of the trace over the fermion structure in the Kubo formula, which 
allows us in this case to freely interchange the Dirac structure of the two vertices. The 
sum of A/22 (a) and A/22 (c) is thus twice the real part of A/22 (a)- Following the procedure of 
Sect. IV A\ we multiply the left hand side of Eq. [51] with u s (k) and get: 

[f,s,h 

^2Re[V^ 2/1 U(A;;p,0 A+(Z)A+(p) 



D F {k) = -P F {k) + 



x 



D%{p) 



(52) 



where D F {k) = u s (k)V F (k). The expression [V^V / 1 *_ +2At ] can be viewed as 1 — ► 2 matrix 
elements squared, where one of the element is an effective amplitude that takes into account 
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the LPM effect. This equation is formally equivalent to the linearized Boltzmann equation 



of Arnold, Moore and Yaffe 



23 



241 ]. including the 1 — > 2 collinear scatterings and the LPM 



effect. To make a closer connection to their results, we again define the deviation from 
equilibrium x l {k) = D l (k)/Yk and follow the procedure in Sect. IV Al The final result is: 



-(l-n F (k°))I F (k) 



d 3 p d 3 l 
(2tt) 3 (2tt) 3 



f,s,h 



>2 fiivm 



(k;p,l) 



x [n F (k°)(l ± n v (p°))(l ± n m (l°)) (*'(*) - xl(p) ~ xL(l))} 



d 3 p d 3 l 



f,s,h 



(2tt)V 4) (A; + p-Z) ^2Re[V^y;_ 



t {k,p;l) 



(27r) 3 (2vr) 3 v ~' v " v '" ' * ^—~L'*"i-*fJ™v 
x [n F (k°)n v (p°)(l ± n m (l°)) (**(*) + X M - xln(l))} (53) 



This last equation is identical to the linearized Bo 



menta non-integrated of Arnold, Moore and Yaffe 23j, |24j] obtained using kinetic theory. 



tzann equation with the transverse mo 



VI. CONCLUSION 



In this paper, we have derived the integral equations needed for the calculation of electrical 
conductivity in QED starting from basic quantum field theory. A visual summary of our 
calculation is presented in Fig. [TTJ With power counting arguments and a constraint on the 
ladder kernel due to the Ward identity, we have included all the necessary rungs for a leading 
order result. Our calculation includes rungs corresponding to 2 — > 2 and 1 — >• 2 collinear 
scatterings, including the LPM effect. An important point of our method is that since we 
used the Ward identity from the onset, gauge invariance is explicitly enforced throughout 
the calculation. Specifically, the constraint on the ladder kernel tells us which rungs must 
be included in the integral equation for each self-energies resummed in the propagators. We 
have finally shown the equivalence between our results and the ones of Arnold, Moore and 
Yaffe obtained from an effective kinetic theory. 

As a final remark let us mention that, even if the present paper is restricted to the 
case of electrical conductivity in hot QED, the Ward identity constraint can in principle be 
applied to other gauge theories and other transport coefficients. Of course, to generalize our 
calculation to non-abelian gauge theories is challenging, because of the presence of ghosts and 
of an even subtler power counting. This is the subject of future research. The shear viscosity 
case is also more complicated than the electrical conductivity case, due to the presence of 
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both fermionic and bosonic effective vertices leading to coupled integral equations. The 
extension of our calculation to shear viscosity is in progress 38] . 
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FIG. 16: Momentum convention used for a) The summation formula (|A2p and the integral equa- 
tion (|Alip b) The general 4-point function (jA8j) 



APPENDIX A: ANALYTIC CONTINUATION OF THE INTEGRAL EQUATION 



In this appendix we show in more detail the steps needed to go from the integral equa- 
tion in Euclidean space (1231) to the one in Minkowski space (124"|) . As is well known from 
finite temperature field theory (see for example 50J), one must do the sum over Matsubara 
frequencies before doing the analytic continuation from imaginary energies to real energies. 
An elegant way of doing this sum is presented in Ref. 18| . Briefly, the method gives a 



formula for the summation of expressions of the type 



(ii>k, iv p , iv q )GF{iVp)V F (iv p , iv v — ih> q )GF{iv p 



(Al) 



where Gf(^p) and Gp{iv v — ii/ q ) are resummed propagators (the momentum convention is 
shown in Fig. [Pol) and the iu's are discrete frequencies. This kind of expression is precisely 
what appears in the integral equation (|23]h since the two propagators on the "side rails" are 
resummed (to regularize pinch singularities). The idea is to replace the sum over discrete 
frequencies by a contour integration with a contour encircling those frequencies; the contour 
is then deformed to encompass the branch cuts coming from the resummed propagators and 
other possible poles contained in F(z). The resulting summation formula is 18j 



T F K) = T Yl n B / F (zi)Res(F, 

even/odd u p poles 



Z = Zi 



(2m) 



it-b/f (0 Disc (F) (A2) 



where £ is a parameter specifying the position along the branch cut. 

To be able to do the Matsubara sum present in Eq. (|23j) . we need to know the analytic 



structure of Vp and M.. As argued in 



18] from induction, T>p has only the singularities 
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of the product of Gf(w p )Gf(w p — iv q ) and thus does not contribute to the sum. On the 
other hand, the M. factor contains poles and must be considered in the sum. The problem 
is that we do not know a priori the precise form of the kernel A4, since this is the object 
we want to put a constraint on. The way around this difficulty is to replace M. with the 
spectral representation of a general 4-point function, of which we know the pole structure 
without specifying its exact form. The program is thus to express the kernel Ai in spectral 
form, insert it in the integral equation ff23l) . do the sum over Matsubara frequencies using 
Eq. (1A2j) and do the analytic continuation towards real enegies. 

In the imaginary-time formalism of finite temperature field theory, the (Euclidean) 4- 
point function is the thermal average of a time-ordered product of four bosonic/fermionic 
fields (in the following, we suppress spatial variables for simplicity): 

M(Ti,T 2 ,T 3 ,n) = (T(0 ai (r 1 )0 a2 (r2)0 a3 (r3)0 a4 (r4))) (A3) 

where the Tj's are imaginary times and the <p ai are bosonic/fermionic fields. The Oj's are 
labels that denote the temporal order of the fields. Since there are four fields, the T-product 
can be expanded into 4! possible time orderings. For simplicity, let us consider a particular 
time ordering {r\ > r 4 > r 2 > r 3 ): 

M aia4a2 a 3 = 0(Ti - T 4 )6(t 4 - T 2 )6>(t 2 - T 3 ) (0 ol (n)0a 4 (^)0a 2 (T 2 )0a 3 (^)) (A4) 

Extracting the time dependence of the Heisenberg operators and averaging with respect to 
an equilibrium density matrix p = e _/3H , we get 

M ai a 4 a 2 a 3 = 0{t X - T 4 )6»(t 4 - T 2 )6»(t 2 - T 3 ) 

x Tr (e-^e*> ai (0)e 6 ^ Va 4 (0)e^" T 4 Va 2 (0)e 6 ^~ T 2 Ua 3 (0)e"^) 

(A5) 

Inserting complete sets of energy eigenstates (with eigenvalues E mi ) between each 
bosonic/fermionic operator and rearranging, we obtain 

M aiaia2a3 = 0(n - r 4 )6»(r 4 - r 2 )6»(r 2 - r 3 ) 

X ^ e -Em 1 (P+T3-T 1 ) e -E m4 {Tl-T4) e -E m2 (T4-T 2 ) e -E rn , 3 (T2-T 3 )^ 

mi 

x (mi|0 Ol (0) |m 4 ) (m 4 |0 a4 (O) \m 2 ) (m 2 \(f>a 2 (0) |m 3 ) (m 3 |0 a3 (0) |mi) (A6) 
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where % runs from 1 to 4. To get other time orderings, we only need to make permutations 
of the di's and the Tj's. Since the integral equation fl23|) is in momentum space, we need 
the spectral representation Ai in momentum space. We thus take the Fourier transform of 

Eq. sm : 



M ai 



141213 



dn / dr 2 I dr 3 I dne luJ1T1 e luJ2T2 e luJ3T3 e lW4T; 



Jo Jo 

11141213 ( 



xMa,^a-,a 3 {n, T 2 , T 3 , T 4 ) 



(A7) 



Note that the Fourier transform goes from to (3 and not — 00 to 00, since the presence of the 
medium breaks Lorentz invariance and imposes periodic (anti-periodic) boundary conditions 
on the imaginary time. Consequently, the frequencies iuji are discrete and are even (odd) for 
bosons (fermions). Doing the Fourier transform and manipulating the expression, we get 



M aiaia2a:j (iui, iu 2 , iu 3 , iuj 4 ) 



<W 2 ,^4 | II / ?S (27r)5(fci + h + h + h) ) x 



(27T) 



71423 



_ [iuj A + iu 2 + iu 3 - 


- (fei + k 2 + k 3 )] [iuj 2 + - 

74231 


- (fe + k 3 )][iu} 3 - 


-k 3 ] 


+ T . ; : — 

[lUJ 2 + IU0 3 + lUJi 


- (k 2 + k 3 + &i)][2Co>3 + 

72314 


- (h + h)][iui 


-h] 


+ T . ; ; 

[lUJ 3 + %UJ\ + ZOO 4 
+ t: : : 


- (A; 3 + ki + fct)][^i + 2Cl> 4 

73142 


- (ki + k 4 )}[iuj A 


-h] 



iuji + ico> 4 + iuo 2 - (£4 + k 4 + fc 2 )][zo;4 + ioo 2 - (/c 4 + k 2 )][iu 2 - k 2 ] 



(A8) 



where 5, 



Ml .M1,M,M< 



functions [59] 



is a Kronecker delta and we have introduced the 4-point thermal Wightman 



satisfying the (generalized) KMS condition 

71234(^1,^2,^3,^4) = (-I)e /3fel 7234l(^l,fc2,^3,fc4) 



(A9) 



(A10) 



Equation ( 1A8I) is identical to the one obtained in 59J and gives the pole structure of one 
possible time ordering and its circular permutations. As before, we obtain other time or- 
derings by permutation. Note that we could have expressed Eq. f]A8|) in terms of 4-point 
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spectral densities (i.e. linear combinations of 4-point thermal Wightman functions) 59], 
but we prefer to leave it in this form. 

With the spectral form (1A8j) (and all the other time orderings), it is possible to do the 
Matsubara sum in the integral equation f[2"3"j) . The summation should be done for all the 24 
time orderings of the 4-point function, but for simplicity we only show the summation of the 
71423 term. Inserting the first term of Eq. (1A8I) in Eq. (|23"1) (with the momentum convention 
shown in Fig. [I6|) . we get 

d 3 p ' 



(2*0 



n 



dlj 

W) 



(2n)5(h + h + h + k) 71423 



X 



G F (iv p ) [Q^V F (iu p , iv p - iu q )] G F (iv v - iv q 



-n'k - l&aWivp + ivh - iv q + hsWiVp + h] J 
where we dropped the inhomogeneous term (which vanishes in the — > limit) and we used 
the shorthand notation Uj = k + lj. Note that the notation is symbolic and does not respect 
the Dirac structure of the equation. Using the summation formula (1A2I) and remembering 
that iv k is odd, iv v is odd and %v q is even in tt, we obtain 



Q^V F (iu k ,iu k 



%v„ 



d 3 p 



n 



dlj 

W) 



(2n)6(h + l 2 + h + k) 71423 



X 



n F (-iv k + iv„ - l 23 )D(-iu k + %v q - l 23 , -iv k - I 



23 1 



-iv k - I 



123 



(2m 



(2m 



■n F (0 



-iuk + iv q - k] 

D(£ + ie,£ + ie 



+ 



n F (-k)D(-l a 



-%v„ 



[-iu k - Z423] - iv q + l 2 ] 



■n F (0 



-iu k - Z423] [f + ie + iu k - iv q + Z23] [f + «e + Z 3 ] 

g(£ ~ jej - ie - ii/g) 

[—iv h - Z423] {£,~ie + iv k - iv q + / 23 ] [f - ie + l 3 ) 

D(j + ie + iv q ,i + ie) 

-iv k - Z 423 ] [f + ie + iz/fc + Z23] [f + it + w q + h] 
D(£ -ie + iv q ,i - ie) 



(A12) 

-iv k - / 423 ] [£ - ie + iv k + Z23] [£ - ie + iu q + l 3 ] ' 
where we used the shorthand notation D(a,b) = G F (a)[Q IM V F (a,b)]G F (b). In the above 
expression, the propagator G F {h) has no ie prescription and is thus undefined. This problem 
can be cured by assigning a small imaginary part to the integrated momenta k, provided 
that it does not violate the delta function constraint. This is allowed since the 4-point 
function is an analytic function and it will remain so no matter how the imaginary parts 



are assigned 



60] . In the present case, two possible choices are Z 3 — ^ Z 3 ± 2ie\ the remaining 
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momenta can be adjusted in many ways so that the delta function constraint is satisfied. It 
is important to note here that this internal assignment of imaginary parts must always be 
"less important" than the ze's coming from the analytic continuation of iu^ and iv q , because 
the latter encode the physics of transport coefficients. Doing the analytic continuation 



k° + 2ie, 



q° + 4ie and doing the integration, the integral equation for Ai 



J 12:! 



becomes: 



lim^P^fc + ie,k" - 



it 



q^O 



d 3 p 



(27T) 



n 



{2ir)8{h + l 2 + h + h) ) 71423 



X 



M-Z23) + n F (-k° - h3))D(-k° - l 23 + ie, —k° - l 23 - ie) 



(A13) 



[k° + U 23 + ie][k° + l 2 -ie] 
where non-pinching pieces have been dropped and the q — > limit has been taken. Note that 
the result is independent of the ie prescription for l 3 . It can be shown by explicit calculations 
that this independence on the assignment of imaginary parts on internal momenta holds for 
all time orderings. 

We can do the same exercise with other time orderings. Collecting all terms together and 
rearranging, the complete Matsubara summed integral equation is (spatial variables are still 
suppressed here and the q — > is implicit) 

d 3 p f dp 



q,V^ 



X 



(2tt) 3 
(1 + e 



;i - n F (p ))Gf(p°)q,V^p°)G^(p^ 



(2*) 



+ (l + e 



-/3k° 



( 72341 - 72314 + 73241 ~ 73214 \ 

(2tt) J (2tt) V [k° + h~ ie] [k° - h + ie] J 

dli f dl2_ ( 71324 + 71342 ~ 73124 ~ 73142 \ 



(2tt) J (2vr) V [k° + h - ie][k° - h + ie] )\ (M4) 
The arguments of the 7's are Tajoaasa* (h, h,P° — h + k°) on the second line and 
701020304(^1? h,P° — h — k°) on the third line. To make sense of this equation, let us compare 
it the the integral equation obtained by Jeon [l3[ (suitably generalized to fermions): 



q^(k°) 



d A p 

WT 4 



n F ( P °)(l + e^°) (M33 + e-^°M 23 ) (1 + e"^)- 
l-n F (p»))GT(p)q,V»(p)Gt(p) 



r, 



X 



l + e /3fc °)n F (p u )(l + e-^ u )M ; 



(3p°\ 



33 



+a + e 



^°)n F (p )(l + e-^)M 23 



(A15) 
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where M33 and M23 represent (completely cut) kernels. Comparing Eqs. (1A14j) and (1A15I) . 
we clearly see that they have the same structure. This is normal, since they both represent 
the same physics. We can thus infer from this that the double integrals over l\ and I2 
in Eq. (1A14I) are the general spectral representation of (completely cut) kernels. The fact 
that kernels (4-point functions) are represented by double integrals (3-point functions) in 
Eq. (I A 141) is normal, because there are only two independent momenta in the q — > limit. 
The similarity between Eqs. (IA14I) and flA15j) also shows that the Euclidean integral equation 
keeps its form when going to Minkowski space. 
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